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We calculate the number density, energy density, transverse pressure, longitudinal pressure, and 
magnetization of an ensemble of spin one-half particles in the presence of a homogenous background 
magnetic field. The magnetic field direction breaks spherical symmetry causing the pressure trans- 
verse to the magnetic field direction to be different than the pressure parallel to it. We present 
explicit formulae appropriate at zero and finite temperature for both charged and uncharged par- 
ticles including the effect of the anomalous magnetic moment. We demonstrate that the resulting 
expressions satisfy the canonical relations, f2 = —P\\ and P± — P« — MB, with M — —dQ./dB being 
the magnetization of the system. We numerically calculate the resulting pressure anisotropy for a 
gas of protons and a gas of neutrons and demonstrate that the inclusion of the anomalous magnetic 
increases the level of pressure anisotropy. 

PACS numbers: 71.10.Ca, 95.30.Tg, 97.10.Ld 



I. INTRODUCTION 

The determination of the bulk properties of a Fermi 
gas in the presence of a magnetic field is important for 
understanding neutron stars and the early-time dynam- 
ics of the quark gluon plasma created in relativistic heavy 
ion collisions. In the presence of a uniform magnetic field, 
both the matter and the field contributions to the space- 
like components of the energy-momentum tensor become 
anisotropic. The degree of pressure anisotropy increases 
as the magnitude of the magnetic field increases. In this 
paper we revisit the calculation of the bulk properties 
of a Fermi gas of spin one-half particles in a uniform 
magnetic field with the goal of unambiguously determin- 
ing the pressure anisotropy from first principles including 
the effect of the anomalous magnetic moment. 

As mentioned above, there is currently considerable in- 
terest in the behavior of matter in the presence of high 
magnetic fields. Neutron stars, for example, are known 
to possess high magnetic fields. More specifically, magne- 
tars [IHI] are believed to have surface magnetic fields as 
strong as 10 14 — 10 15 Gauss. Based on such surface mag- 
netic fields, one could expect magnetic fields in the inte- 
rior of magnetars to be on the order of 10 16 — 10 19 Gauss. 
There have been many previous studies of the effect of 
magnetic fields on neutron stars and magnetars focusing 
on the effect of magnetic fields on the equation of state of 
the matter composing the star including hadronic mat- 
ter, quark matter, and hybrid stars composed of hadronic 
matter with a quark matter core [81448] . 

Among these references some authors have simply as- 
sumed that the system continues to be describable in 
terms of an energy density and an isotropic pressure 
derivable from standard thermodynamic relations, while 
other authors have included the fact that the background 
magnetic field breaks the spherical symmetry of the sys- 
tem. The breaking of the spherical symmetry has two 
distinct contributions: (i) the matter contribution to the 



energy-momentum tensor and (ii) the field contribution 
to the energy-momentum tensor. For charged particles 
the presence of a magnetic field causes the pressure trans- 
verse to and longitudinal to the local magnetic field direc- 
tion to be different, with the level of pressure anisotropy 
increasing monotonically with the magnitude of the mag- 
netic field. On the other hand, for uncharged particles 
one finds that the pressure is isotropic. 

There have been dynamical models of neutron stars 
which have attempted to include the effect of high mag- 
netic fields on the three-dimensional structure of neutron 
stars [49 54 . Some of these studies have self-consistcntly 
included modifications of the general relativistic metric 
necessary to describe the breaking of spherical symmetry 
by the neutron star's magnetic field. However, to the best 
of our knowledge there has not been a study which has 
simultaneously included the general relativity aspects, ef- 
fects of magnetic fields on the equation of state, and ef- 
fects of pressure anisotropy on the static and dynamical 
properties of a high-magnetic-field neutron star. In order 
to complete this program it is necessary to first under- 
stand all sources of pressure anisotropy due to magnetic 
fields. 

Another area in which there has been a significant 
amount of attention focused on the behavior of matter 
subject to high magnetic fields is the consideration of the 
first fm/c after the collision of two high-Z ions in a rela- 
tivistic heavy ion collision. Because of the large number 
of protons in the colliding nuclei, magnetic fields on the 
order of 10 18 — 10 19 Gauss are expected to be generated at 
early times after the initial nuclear impact [55-59 . The 
existence of such high magnetic fields prompted many re- 
search groups to study how the finite temperature decon- 
finement and chiral phase transitions are affected by the 
presence of a background magnetic field. These studies 
have included direct numerical investigations using lat- 
tice quantum chromodynamics (QCD) [60-63] and theo- 
retical investigations using a variety of methods includ- 
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ing, for example, perturbative QCD studies, model stud- 
ies, and string-theory inspired anti-de Sitter/conformal 
field theory (AdS/CFT) correspondence studies [53H57]. 

In order to have more a comprehensive understanding 
of the behavior of matter in a background magnetic field, 
we begin with the basics and study Fermi gases consisting 
of charged and uncharged spin one-half particles includ- 
ing the effect of the anomalous magnetic moment. Many 
of the results obtained here are already available in the 
literature; however, the results for the transverse pres- 
sure including the effect of the anomalous magnetic mo- 
ment have not appeared previously. For sake of complete- 
ness, we present the results for all of the components of 
the matter contribution to the energy-momentum tensor 
with and without anomalous magnetic moment as a point 
of reference for future applications. In this paper we con- 
sider systems at both zero and finite temperature. For 
zero temperature systems, we demonstrate by explicit 
calculation that the grand potential Q = e — /in = —P\\ 
where e is the energy density, n is the number density, 
Pv, is the pressure along the direction of the background 
magnetic field, and [i is the chemical potential. For finite 
temperature systems one also finds that O = — P||. 

We then show that the resulting expressions satisfy 
the canonical relation P± = Pu — MB, where P± is the 
pressure transverse to the magnetic field direction and 
M = —dil/dB is the magnetization of the system. Eval- 
uating the resulting expressions numerically, we demon- 
strate that the magnitude of the pressure anisotropy is 
larger when one takes into account the anomalous mag- 
netic moment, however, as the temperature of the system 
increases the pressure anisotropy decreases. 

The structure of the paper is as follows. In Sec. |H| we 
introduce the basic formulae necessary to calculate the 
bulk properties of an ensemble of particles using quantum 



field theory. In Sec. Ill we present the resulting formulae 
for charged particles with and without anomalous mag- 
netic moment. In Sec. |IV| we present the corresponding 
formulae for uncharged particles. In Sec. [V] we compare 
the numerical evaluation of the transverse and longitu- 
dinal pressures. In Sec. |VI| we present our conclusions 
and an outlook for the future. Finally, in Apps. |A"|and|B| 
we present a quantum field theory derivation of the nec- 
essary components of the energy-momentum tensor for 
charged and uncharged particles. 



II. GENERALITIES 

In the presence of fields, the energy-momentum tensor 
can be decomposed into matter and field contributions 
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If there is only a background magnetic field B point- 
ing along the z-direction, then the field contribution to 
the energy-momentum tensor takes the form T^ ds = 



diag( J B 2 /2,B 2 /2,S 2 /2,-B 2 /2)Q Since this contribu- 
tion is well-understood, we do not spend more time dis- 
cussing it in this paper. Instead, we focus on T^ tter for a 
system composed of spin one-half fcrmions subject to an 
external magnetic field. In what follows, the bulk prop- 
erties of the system (energy density, transverse pressure, 
etc.) are understood to specify the components of Tetter 
in the local rest frame of the system. 

The matter contribution to the bulk properties of a sys- 
tem can be expressed in terms of the one-particle distri- 
bution function /. We consider a single particle type with 
mass m and charge q and sum over the spin polarizations. 
The results obtained can be straightforwardly extended 
to a system consisting of multiple particle types. We 
present a derivation of the necessary components of the 
energy-momentum tensor in Apps. |A"|and|Bl Summariz- 
ing the results, one finds that the local rest frame number 
density, energy density, longitudinal pressure, and trans- 
verse pressure can be expressed in terms of the following 
integrals of the one-particle distribution function 
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where we have singled out the z (parallel) direction for 
future application, fr?{y) = (a/to 2 + fc 2 — skB) 2 , kj_ is 
the (discretized) transverse momentum, ^ s represents 
a sum over spin polarizations, k represents the anoma- 
lous magnetic moment, and J, is a properly normalized 
(sum-)integration over momenta which we will define sep- 
arately for charged and uncharged particles. For charged 
particles with vanishing anomalous magnetic moment, 
the expressions above were first derived in Ref. [8]. For 
charged particles with finite anomalous magnetic mo- 
ment the expressions for the number density and energy 
density above were first derived in Ref. llj. Here we 
extend the treatment to include uncharged particles and 
independently compute the transverse and longitudinal 
pressures in the case of finite anomalous magnetic mo- 
ment. 

We note that in order to include interactions, one 
should use the interaction-corrected expression for the 



1 This is the form in Heavisidc-Lorentz natural units. In Gaussian 
natural units, when converting the magnetic field to GeV 2 , the 
magnetic field is increased by a factor of vJtt and the components 
of the energy momentum-tensor are divided by a factor of 4ir to 
compensate, e.g. £g = B 2 /8tt. 
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particle's dispersion relation. In the mean-field approxi- 
mation, this amounts to including corrections to the bare 
mass of the particle being considered, e.g. m — > m* . The 
resulting effective mass can depend on the chemical po- 
tential and temperature. In what follows we indicate the 
effective mass of the particle as m assuming that inter- 
action corrections could be absorbed into the mass0 



III. CHARGED PARTICLES 

In the presence of a uniform external magnetic field 
pointing in the z-direction, the transverse momenta of 
particles with an electric charge q are restricted to dis- 
crete Landau levels with k\ — 2v\q\B where v > is an 
integer [S5] and one has 
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where the sum over n represents a sum over the dis- 
cretized orbital angular momentum of the particle in the 
transverse plane. For spin one-half particles the orbital 
angular momentum n is related to v via 



v = n + 



1 
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where s = ±1 is the spin projection of the particle along 
the direction of the magnetic field and q is the charge |j 
An additional consequence of the quantization is that 
the total energy of a charged particle becomes quantized 



E=\jk 2 z + ((m 2 + 2v\q\B) 1 / 2 - snB) 2 . 
— \Jk 2 z + fh 2 {v) , 



(8) 



where k = K^/ijv with «j being the coupling strength 
for the anomalous magnetic moment times the magneton 
and fh 2 {v) = (^/m 2 + 2v\q\B — skB) 2 . 



A. Zero temperature 

At zero temperature the one-particle distribution func- 
tion is given by a Heaviside theta function 

f(E)=Q(»-E), (9) 

where fi is the chemical potential (Fermi energy). 



2 In the following, spherical symmetry is broken by a uniform mag- 
netic field. Due to this, the effective mass could, in principle, also 
depend on the angle of particle momentum relative to the mag- 
netic field direction. We do not take this possibility into account 
in this work. 

3 The present calculation is valid only for spin one-half particles. 
Spin zero, one and three-half particles, described respectively 
by the Klein-Gordon, Proca and Rarita-Schwinger equations are 
affected differently by the magnetic field | 86l I89| . 



1. Zero anomalous magnetic moment 

We begin by considering the case with no anomalous 
magnetic moment, i.e. K = 0. In terms of the chemical 
potential, //, the maximum k z is defined via @ 



kzA") = VM 2 -2%|S-m 2 . (10) 
In addition, in the sum over the Landau levels one must 



guarantee that the quantity under the square root in ( 10 ) 
is positive. This requires fh 2 < /i 2 which results in 



V < Vr, 



2\q\B 



(11) 



where [x\ — max{n G Z | n < x} is the largest integer 
less than or equal to x. 

Using the above, we can write down an expression for 
the number density using ([2| and ^ to obtain 



V<Vrr 



E 



\q\b 
\q\B v 

2tt 2 ^ 
2tt 2 ^ 

s=±l n=0 



dk z Q(fi-E), 



dk z 



(12) 



Note that the upper limit on the n sum is set in terms 
of the maximum Landau level and that v depends on n 
and s via Eq. Note that the k — degeneracy factor 
for a given Landau level is automatically taken into ac- 
count by the dual sum over spin and angular momentum. 

Similarly, one can evaluate the energy density to obtain 
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Next, we consider the parallel pressure P\\ and obtain 

p \i\b v f k » F „ k 2 
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Note that using (12 1, (13), and (14 1 it is straightforward 
to see that e + P» = [in and hence = e — fin = — Pii . 
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Finally, we consider the transverse pressure Pj_ and 
obtain 
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Numerically the results for Py and Pj_ are different for 
any value of P; however, they only become significantly 



different for very large P. Using Eq. (Ill, for example, 
we see that when B > (fi 2 — m 2 )/2|g[7~only the lowest 
Landau level contributes to the sums and one obtains 
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where kp = \J — m 2 . The transverse pressure on the 
other hand vanishes in this limit 



lim Pi = . 
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In the case that there is no anomalous magnetic mo- 
ment, a relationship between Pm and Pj_ can be es- 
tablished by evaluating the magnetization of the sys- 
tem M ee -dfl/dB = dP\\/dB $U\. Performing the 
necessary derivative^] of the parallel pressure one finds 
M = (P|| - PjJ/P. Rearranging gives Pj_ = P\\ - MB 
which is the canonical relationship one finds in the liter- 
ature between the transverse and longitudinal pressures. 



2. Nonzero anomalous magnetic moment 

We now turn to the case of nonzero anomalous mag- 
netic moment. In this case the expressions for k Zt p and 
c max must be adjusted to 



k z ,F = vV - (( m2 + 2H<?|£) 1/2 - skB) 2 , (18) 
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The transverse pressure, however, is modified when 
there is a non- vanishing anomalous magnetic moment 
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Evaluating the magnetization one obtains in this case 
OB ~ B ' 2tt 2 ^ 
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So one finds once again P± = Pm — MB. 



B. Finite temperature 

We now turn our attention to the case of a finite tem- 
perature ensemble of charged particles. In this case the 
distribution function is 



f±(E,T,fi) 



1 
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where /+ describes particles, /_ describes anti-particles, 
and [i is the chemical potential. 



1. Zero anomalous magnetic moment 
We begin with the number density 

MP 
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With these two modifications Eqs. (12), (13), and ([14) ^ E = f k 2 + rh 2 {v) with 



are unchanged, but one should note that ^ max now de- 
pends on the spin alignment s. 



2v\q\B. Introducing the variable x — E^f^i we can rewrite 
k z = yj(x ± fi) 2 — fh 2 (y) and using dk z = (x ± fi)((x ± 
^) 2 — TO 2 (^)) _1 / 2 dx one obtains 



4 Formally one should use left or right derivatives in the vicinity of 
magnetic field magnitudes where fmax changes under infinitesi- 
mal variation. 
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Next, we consider the energy density. Using the same 
change of variables as before, one obtains 



\q\B 
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Similarly, one obtains for the longitudinal pressure 



P \\,± = 
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dx ^{x±fi) 2 -fh 2 (v)f± (x, T, 0) . 

')=FM 

(26) 

Finally, one obtains for the transverse pressure 



P±,± = 
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(27) 

Next we consider the magnetization obtained from 
M = dP\\/dB. In order to do this we apply the fun- 
damental theorem of calculus 
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Using this we can evaluate the derivative of the integral 
appearing on the second line of ( 26 ) 
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dx yj{x±p) 2 -m 2 {v)f ± (x, T, 0) = 
where we have used the fact that in the case at hand the 
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first term on the right-hand side of (28) is zero. Using 
fhdrh/dB = ^dm 2 /dB = \q\u we can obtain finally 



M+ = 



dP\\ 



P\\ 



P 



-L,± 



dB 



B 



B ' 



(30) 



which is the canonical relation between the transverse 
pressure, the longitudinal pressure, and the magnetiza- 
tion. Rearranging we obtain P±.± = P\\.± — M±B be- 
tween the perpendicular and parallel pressures at finite 
temperature in the case that there is no anomalous mag- 
netic moment. 



2. Nonzero anomalous magnetic moment 

As was the case at zero temperature, when including 
the anomalous magnetic moment, the primary thing that 
changes is the mass fh 2 (v) — {^Jm 2 + 2v\q\B — skB) 2 . 



With this change, the expressions for n±, e±, and P\\,± 
given in Eqs. (24), (25), and (26), respectively, are un- 
changed. For the transverse pressure, however, one must 
include additional terms 
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In addition, when including the anomalous magnetic mo- 
ment, the magnetization has a different form since 
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Once again we see that P±± — Pu ± — M±B. 



IV. UNCHARGED PARTICLES 

In the case that the particle being considered is un- 
charged, then one does not obtain discrete Landau levels 
and, as a result, 



d 3 k 
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in Eqs. ([2])-([5|. Prior to proceeding with the calculations, 
we note that for uncharged particles one has 



m 2 + fc 2 ^ — skB 



(35) 



A. Finite temperature 

We first consider the general case of uncharged par- 
ticles at finite temperature including the effect of the 
anomalous magnetic moment. The derivation necessary 
is performed in App. [B] Here we summarize the re- 
sults and list the contributions from particles and anti- 
particles. The resulting expression for the number den- 
sity is 
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where k = \/E 2 — (m — skB) 2 . The energy density is 
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The longitudinal pressure is 
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The transverse pressure is 
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Finally, we obtain the magnetization 
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We see that the magnetization vanishes when k — > 0. In 
addition, with these expressions one finds P±,± = P\\,± — 
M±B. 



where to = to — skB, kp — \J li 2 — to 2 . Similarly the 
energy density can be obtained in this limit 
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And the longitudinal pressure can also be easily obtained 
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Using the derived expressions for n, e, and Pm one can 
show that e + P|| = /m is satisfied explicitly. The trans- 
verse pressure is 



S=±l 



(2fi 2 - 5to 2 - UsnBrh - \2{snBf 



-3m 2 (m + 2 S KB) 2 \og l kF + ^ 



• (44) 



Finally, evaluating dP^/dB one obtains the magnetiza- 
tion in this case 
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From this result we can once again verify that P± 
Pn - MB. 



B. Zero temperature 

In the zero temperature limit there is only a particle 
contribution since liniT->o f-(E, T, fi) — for E > and 
lim T ^ f+{Ei Pi m) = QCA 4 — £■)• Using the results listed 
in the previous subsection one finds for the number den- 
sity 
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V. NUMERICAL RESULTS 

In this section we present numerical evaluation of the 
transverse and longitudinal pressures derived in the pre- 
vious section. For the numerics that follow we will as- 
sume (i) a gas of protons with a mass to = m p — 
0.939 GeV, electric charge q — +e, and an anomalous 
magnetic moment of K = rtp/ijv = 1-79 • e/{2m 2 ) — 
0.288633 GeV" 1 in Heaviside-Lorentz natural unitsnand 



5 In Gaussian natural units one has fJ.M = 0.0454871 GeV 1 which 
is the Heaviside-Lorentz value divided by \/4tt. Note that if one 
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n [fm" 3 ] 

FIG. 1. (Color online) Transverse and longitudinal pressures 
of a zero temperature gas of protons as a function of the num- 
ber density. Results include the effect of the proton anoma- 
lous magnetic momentum. 




-1 



0.2 0.4 0.6 0.8 1 
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FIG. 2. (Color online) Ratio of transverse and longitudinal 
pressures of a zero temperature gas of protons as a function 
of the number density. Results are shown with and without 
the effect of the proton anomalous magnetic moment. 

(ii) a gas of neutrons with a mass m = m n = 0.939 GeV, 
electric charge q — 0, and an anomalous magnetic mo- 
mentum of k — K n HN — —1-91 • e/(2m„) = —0.307983 
GeV -1 [3S]. In all cases shown we consider a magnetic 
field magnitude of 5 x 10 18 Gauss. 

In Fig. [I] we plot the transverse and longitudinal pres- 
sures of a zero temperature gas of protons including the 



uses Gaussian natural units, the magnetic field in GeV 2 is scaled 
by a factor of v47r compared to the corresponding Heaviside- 
Lorentz magnetic field. As a result the product of /^jvB is inde- 
pendent of the convention chosen. 
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FIG. 3. (Color online) Magnetization of a zero temperature 
gas of protons times the background magnetic field. Result 
includes the effect of the proton anomalous magnetic moment. 




FIG. 4. (Color online) Ratio of transverse to longitudinal 
pressure of a gas of protons as a function of the net proton den- 
sity for four different temperature T = {0, 10, 30, 500} MeV. 
Results include the effect of the proton anomalous magnetic 
moment. 



effect of the anomalous magnetic moment. The cusps in 
the curves correspond to threshold crossings for the max- 
imum Landau level. As can be seen from this figure, the 
transverse and longitudinal pressures are not equal. In 
addition, one can see from the figure that at low densi- 
ties the transverse pressure is negative at low densities 
when there is a non-vanishing anomalous magnetic mo- 
ment, while the longitudinal pressure remains positive at 
all densities. 

In Fig. [2] we show the ratio of the transverse to longi- 
tudinal pressures for a zero temperature gas of protons 
with and without the effect of the anomalous magnetic 
moment. In both cases we once again see cusps indicative 
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FIG. 5. (Color online) Ratio of transverse and longitudinal 
pressures of zero temperature gas of neutrons as a function 
of number density. Results are shown with and without the 
effect of the neutron anomalous magnetic moment. 
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FIG. 6. (Color online) Ratio of transverse and longitu- 
dinal pressures of a gas of neutrons as a function of the 
net neutron density for four different temperature T = 
{0, 10, 30, 500} MeV. Results include the effect of the neutron 
anomalous magnetic moment. 



of Landau level crossings and a vanishing transverse pres- 
sure at low densities. From this figure we also see that 
including the anomalous magnetic moment enhances the 
pressure anisotropy. 

In Fig. [3] we plot the background magnetic field times 
the magnetization of a zero temperature gas of protons 
obtained via Eq. ( 21 ). We note that there are two distinct 



sets of cusps visible in Fig. [3] This is due to the fact that, 
when the effect of the anomalous magnetic moment is 
included, there are two different Landau level thresholds 
for particles with spins aligned or anti-aligned with the 
background magnetic field. 



In Fig. [4] we plot the ratio of the transverse pressure to 
the longitudinal pressure of a gas of protons as a function 
of the net proton density (particle minus anti-particle) for 
T = {0, 10, 30, 500} MeV. As can be seen from this figure, 
as the temperature is increased, the cusps associated with 
Landau level crossings are diminished and the level of the 
pressure anisotropy also decreases. The highest tempera- 
ture shown T = 500 MeV is on the order of those initially 
generated in relativistic heavy ion collisions at CERN's 
Large Hadron Collider. As we see, at these high tem- 
peratures the pressure anisotropy for charged particles is 
quite small, < 1%. However, it should be noted that as 
the system cools, the pressure anisotropy increases. 

We consider next the case of neutral particles, focusing 
on a specific example of a gas of neutrons. In Fig. [5] we 
plot the ratio of the transverse to longitudinal pressures 
of a gas of neutrons as a function of the neutron den- 
sity with and without the effect of the neutron anoma- 
lous magnetic moment. This figure shows that without 
the anomalous magnetic moment the pressures are com- 
pletely isotropic; however, when there is a non-vanishing 
anomalous magnetic moment the pressure anisotropy can 
be quite sizable. In Fig. [6] we show the ratio of the to- 
tal particle plus anti-particle transverse to longitudinal 
pressures. This figure shows that as the temperature of 
the system increases, the amount of pressure anisotropy, 
again, decreases. 



VI. CONCLUSIONS AND OUTLOOK 

In this paper we have revisited the calculation of the 
matter contribution to the energy-momentum tensor of 
a Fermi gas of spin one-half particles subject to an ex- 
ternal magnetic field. We considered both charged and 
uncharged particles with and without the effect of the 
anomalous magnetic moment. For zero temperature sys- 
tems we demonstrated through explicit calculation that 
the resulting energy density, number density, and longi- 
tudinal pressure satisfy £ + P\\ = ^n. Using the standard 
definition of the grand potential £1 = e — [in allowed us 
to see that, in all cases investigated, the grand potential 
is related to the longitudinal pressure via Q = — Pu in 
agreement with previous findings in the literature. 

We point out that many of the results contained herein 
are known in the literature. However, the results ob- 
tained for the transverse pressure in the case that there 
is a non-zero anomalous magnetic moment are new. In 
addition, we have presented in two appendices an ex- 
plicit derivation of the energy-momentum tensor taking 
into account the anomalous magnetic moment. Using 
the results obtained, we demonstrated that the standard 
relationship, P± = Pu — MB, between the transverse 
pressure, longitudinal pressure, and magnetization of the 
system holds in all cases considered. 

The resulting formulae for the bulk properties can be 
applied to both zero temperature and finite temperature 
systems and hence could be useful in understanding the 
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impact of high magnetic fields on the evolution of proto- 
neutron stars, proto-quark stars, and the matter gener- 
ated in relativistic heavy ion collisions. Applying the de- 
rived formulae to a system of protons we found that there 
can exist a sizable pressure anisotropy in the matter con- 
tribution to the energy-momentum tensor which could 
have a phenomenological impact. Additionally we found 
that as the temperature of the system increases, the pres- 
sure anisotropy decreases. This is primarily due to the 
fact that increasing temperature allows higher Landau 
levels to be partially occupied and hence reduces the dis- 
crete effects one sees at zero temperature. For uncharged 
particles Landau quantization does not play a role and 
instead any pressure anisotropy exhibited comes from a 
non-vanishing anomalous magnetic moment. Once again 
as the temperature increases, the pressure anisotropy is 
reduced. This effect is due to the fact that as the tem- 
perature increases high momentum modes become highly 
occupied which causes momentum terms in the energy to 
dominate over those associated with the anomalous mag- 
netic moment. 

We note that although we presented results applica- 
ble to the case of a single particle type, the resulting 
formulae can be easily applied to the case of a system 
composed of multiple particle types. Since the contribut- 
ing particles may have different pressure anisotropics de- 
pending on the sign and the magnitude of the anomalous 
magnetic moment, one must take care to sum over all 
particle types subject to the necessary conservation laws 
prior to making quantitative statements about the phe- 
nomenological impact of magnetic-field induced pressure 
anisotropies on dense matter [91j . Finally, we emphasize 
that although the numerical results shown in the results 
section assumed a particular magnetic field amplitude, 
the analytic results derived herein are completely gen- 
eral and as such can be applied to assess the impact of 
magnetic fields on the bulk properties of matter in a wide 
variety of situations. 
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Appendix A: Energy-momentum tensor 

In this appendix we derive the energy-momentum ten- 
sor including the effect of the anomalous magnetic mo- 
ment. For this purpose we will use the method of metric 
perturbations which allows one to most efficiently com- 



pute a symmetric and gauge-invariant energy-momentum 
tensor. The starting point is the following relation be- 
tween the variation of the action and the energy momen- 
tum tensor 



SS 



(Al) 



where g = det(g M „). We proceed in the standard way by 
writing the action in terms of the Lagrangian density, 
varying the metric, identifying the energy-momentum 
tensor by comparison with ( Al ), and finally taking g^ v — ¥ 
r^ v where r\^ v — diag(l, — 1, — 1, — 1) is the Minkowski- 
space metric. 

We begin with the curved-space Lagrangian density 
for a spin one-half fermion with charge q in the presence 
of an external magnetic field including the effect of the 
anomalous magnetic moment 

C = i>{ilp -m+ -/wrT^ - \f^F^ , (A2) 



where n is the anomalous magnetic moment and, as 
£> M = i(% - %) + T M + iqA„ with 



usual, ju 



r M being the spin connection which is zero in flat space, 
and <7^ lV — i[7 M , 7"]/2. This allows us to write the co- 
variantized action as iS = J d A x \J—g C = S m + Sf with 



S m — 



rD^g^+g^) 



S 



f 



+ ^K<T al3 F jS (g Q1 + g ia ){gfis + 9sp) 
\ J d 4 X ^F^F^g al g ps , 



i>, (A3) 
(A4) 



where we have split the action into matter and field con- 
tributions and used the the fact that the metric tensor 
is symmetric to explicitly symmetrize the matter contri- 
bution. First, we evaluate SS making use of the identity 
5y/—g = — \\J—gg ilv Sg >liV . Note, importantly, that the 
gamma matrices themselves depend on the metric and 
therefore one needs to take into account their variation 
under metric variation. The variation can be computed 
with [92] or without [93] the use of vierbeins. Comput- 
ing the variation, identifying T >iV , and taking the limit 
qiiv _^ one nnc j s t ne following expressions for the 
matter and field contributions to the energy-momentum 
tensor in flat space 



f w 



f 



K(<7" a n +(T Ua F^ a ) 



i> - <"£ m , (A5) 

Tf = —F m F v a - ifLj , (A6) 

where C m and C / are the matter and field contributions 
to the Lagrangian density corresponding to the first and 
second terms in Eq. ( A2 1 , respectively. 
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Appendix B: Matter contribution to 



= B(5» x 5 vy - 8 UX 8^ V ) and as a result 



In this appendix we derive expressions for the energy- 
momentum tensor in a uniform background magnetic 
field. We focus on the matter contribution since the field 
contribution (A6| is standard. In the rest of this ap- 



pendix we can therefore ignore the pure gauge field term 
in the Lagrangian. In flat space the Lagrangian density 
for a spin one-half fcrmion with charge q in the presence 
of an external magnetic field including the effect of the 
anomalous magnetic moment is 

C = $(0-m+~Ka' a, F MV )i!>, (Bl) 



where n is the anomalous magnetic moment and, as 

— > i — 
= |(4i - dp) + iqA^, and 

jJ,V 



usual, 



a**" = i[7 M , 7"]/2. The equations of motion for ip and 
tp can be determined using 



dC 
dC 



dC 



dC 



which result in 



(i$ -q4-m+ -KO^F^ = , 

+ 4>(q4 + m- -Ktr^Fp,) = . 



(B2) 

(B3) 
(B4) 



We note for application to the calculation of T^ v that 
if we multiply the first equation from the left by ip we 
obtain C = 0. This demonstrates that the matter La- 
grangian density vanishes when evaluated with solutions 
which obey the equations of motion. This allows us to 
simplify Eq. | A5 1 to 



T vv = il> 



{^D v + YD") 



-K(a» a F v a + <T v<x F» a ) 



i>, (B5) 



To evaluate the necessary statistical average of T" v we 
first need to solve the equations of motion in order to de- 
termine the energy eigenvalues and spinors. The spinors 
and energy eigenvalues are available in the literature [Ml - 
197] ; however, we review the derivation for sake of com- 
pleteness and then use the resulting spinors to evaluate 
the statistical averages of T Ml/ . We note that the result- 
ing spinors have been expressed in various different forms 
in the literature. We present a specific compact form for 
the spinors, however, we have explicitly verified that us- 
ing the forms of the spinors presented in Refs. [91H97] 
yields the same final results. 

As in the main body of the text, we choose the mag- 
netic field to point along the z-direction. Choosing 
the vector potential to be A" — (0, —By, 0,0) we have 



^na^F^ = iKB-f-r =kb( a * ^ \ = nBS 3 . (B6) 



Next we write Eq. ( B3 ) in Hamiltonian form by searching 
for static solutions of the form ip = e~ iEt ^(x) which 
results in the Dirac-Pauli equation [9S] 



(a ■ TV + 7°m - kBj°S 3 )^ = 



(B7) 



where a = 7°7 and 7r = — iV — qA. 

Here we are interested in the diagonal components of 
T" v which for a constant magnetic field are given by 

7-° = ^(n £>V' (B8) 

T xx = $ {i 1 x D x - KBa xy ) ip , (B9) 
T vv = $ (i~f v D y - nBa xv ) ip , (B10) 
T zz = i){i"i z D z )4>. (Bll) 

1. Charged particles 

We now search for the solution of the Dirac-Pauli equa- 
tion for charged particles. Based on the structure of the 
equation, we begin by making an ansatz for the bi-spinor 
# of the form *(x) = e ik * x e ik **v$ (y) with [M] 



C24>u-i{y) 

\ c 4 I/ _ 1 (y) 



(B12) 



where 



1 s q 



(B13) 



with n = 0, 1, 2, • • • The constants Cj above implicitly 
depend on the spin alignment s = ±1. The functions 4> n 
are given by 



where the variable £ is 



II - 



qB 



(B14) 



(B15) 



n > is an integer, H n is a Hermite polynomial, and 
N n = (gS) 1 / 4 (Y / 7r2"n!) -1 / 2 is a normalization constant 



6 When v = there could be an issue with the Hermite functions 
with index u — 1 not being well defined; however, as we will 
show below, in this case one finds that the coefficients vanish 
identically. 
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which ensures dy <j>\(y) = 1. Inserting this ansatz 
and simplifying the Dirac-Pauli equation, one obtains 



X = E X , 



(B16) 

where x = ( c i c 2 C3 C4) T and k v = \Jl\q\Bv. Evaluating 
the determinant of the matrix on the left we obtain the 
energy eigenvalues 



/ m 


- kB 





k z 







m + kB 


k v 


-k z 




k z 


k v 


—m + kB 





V 


ku 


~k z 





— m — kB 



E s = ±y/k 2 z + (\-SKB) 2 . 



(B17) 



where A = ^Jm? + k 2 .. The choice of an overall pos- 
itive sign for the energy eigenvalue above corresponds 
to particle states and the negative sign to anti-particle 
states. Without loss of generality we can focus on the 
positive energy states and, in the end, extend the result 
to include the necessary contribution from the negative 
energy states. 

The resulting positive energy eigenvectors are 



X 



(») 




V2Ao~# 



(B18) 



where a s = E s — kB + sX and j3 s = X + sm. The 
overall normalization of the state is fixed by requiring 
that dyu^^x.) ti|j'(x) = 2E s 5 rs 5 nm . The general 
quantum state for positive energy states can now be con- 
structed 



s=±l 



E Z> s (k)« (s) (k)e^ 



(B19) 



where b s (k) is a particle creation operator which obeys 
{6 r (p), 6j(k)} = (27r)S rs 5 nm S( Pz - k z ) , (B20) 
k = (n, k z ) with n = 0, 1, 2, • • • , k = (E k ,k x ,0, k z ), and 



dk z 1 

27 TUT 



(B21) 



Note that the factor of y / 2E k in the denominator above 
is fixed by the spinor normalization used above. 



man 



H = 



i / %ftd t i> 

J X 

*E 



6j(p)u (r)t (p)e^" 

x J p J k 

& a (k)uM(k)(-*£ fc )e- < *'** 



(B22) 



where J = J d 3 x. Using the orthonormality relations 
listed above one finds 



n = — e e r t^^^ (k)6s(k) - (b23) 



\q\B 

2tt 



s=±l n 

We can now compute the thermal average of the energy 
using the density matrix p 

p = e -> iH+aN , (B24) 

where is fi — l/T is the inverse temperature, a = (3p with 
p being the chemical potential, H is the Hamiltonian 
operator, and 



N = 



\q\B 
2tt 



EEri^W' (B25) 



s— ±1 n 



is the number operator. The statistical average of the 
Hamiltonian operator gives the energy density 



(H) = 



Tr[pH] 
Tr[p] 



(B26) 



Using the Baker-Campbell-Hausdorff formula one ob- 
tains 



(6t(k)6 a (k)) = <6 a (k)6t(k))e-^-"), 



(B27) 



which, upon application of the anti-commutation rela- 
tions for the creation operators, gives the Fermi-Dirac 
distribution for particles 



<& f s (k)& s (k)> 



1 



e 0(E k -n) + 1 



f+(E k ,T,p). (B28) 



With this we obtain our final expression for the particle 
contribution to the energy density 



" = "J" E E 

s=±l n 



2tt 



dk z 
2tt 



Ekf+(E k , T, p) ■ 



(B29) 

Note that if one includes the anti-particle states, one 
must normal order the Hamiltonian operator prior to per- 
forming the statistical average. 



a. Energy Density 

To determine the energy density, we begin by evaluat- 
ing the 00 component of the energy-momentum density 
which is equivalent to the Hamiltonian density T no = 
% = iip^d t ip- Integrating over space gives the Hamilto- 



b. Number Density 

Based on the above discussion, the number density can 
easily be seen to be given by 



n=(N) = 



\q\B 

2tt 



/c 
_ . -< 



s— ±1 n 



dk z 
2tt 



f+(E k ,T,p). (B30) 
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c. Longitudinal Pressure 

We now consider the longitudinal pressure which is 
given by P|| = (J X T ZZ ) with 



T zz = g [VT*^ - (5 Z ^) 7 Z V] ■ (B31) 
Plugging in the explicit forms for the spinors we have 



x J p J k 



x 7 ° 7 z 



&j;(p) U (r)t (p)e^ a;A ' 

fe s (k) u ( s )(k)(ifc z ) e -^ a 

&J(p)u (r)t (p)Hf>*)e^ X *" 
x 7 V ^(kJuWCkJe-* 6 "*"! }.(B32) 



Evaluating the x and p (sum-)integrals, making use of the 
orthonormality relations and then taking the statistical 
average gives 



l\q\B 



2 2tt jtt^J-oo **E k 



/OO 
d 2 /[ U ( s ) t (k) 7 7 z u( s )(k)], (B33) 
-oo 

where we have used the fact that (6£(k)b s (k)) vanishes 
unless r = s. Next we need to evaluate the spinor con- 
traction 

/oo 
dy u (s)t (k) 7 V " (s) (k) = -2k z , (B34) 
-oo 



which follows from the explicit form of the spinors ob- 
tained previously. Using this and rewriting the statis- 
tical average of the number operator as a Fermi-Dirac 
distribution, one obtains 



Z7F s =±l „ -f-oo ^ 



d. Transverse Pressure 

We finally turn our attention to the transverse pres- 
sure. By rotational symmetry, P± = (f T vv ) = 
(f T xx ). Choosing the former, which is somewhat easier 
to evaluate, we should integrate and statistically average 

jvv = - [^vqv^ _ (d^)^] _ nB^o xv i\). (B36) 



Plugging in the explicit forms for the spinors we have 

&t(p)w (r)t (p)e^" 



J X ^ r s J X J p J K 



x 7 V 



b s {k)d v u^(k)e- ik ^ 
x 7 V k(k)uW(k)e- <fc '* x " 

r s J x Jp Jk 

& s (k)uW(k)e-^ x "l . (B37) 



x~f°o- xy 



Evaluating the x and p (sum-)integrals making use of the 
orthonormality relations and then taking the statistical 
average gives 



dy u(^ t (k) 7 °7 ? 'a s/ u (s) (k) 

-a^ (s)t (k) 7 V« (s) (k) 

-kB j dyu^ikj^a^u^ik)^. (B38) 

Integrating by parts one finds that the second term con- 
tributes the same as the first. Using the explicit repre- 
sentation of the spinors obtained above one finds 

^^(khV^u^k) 

i 

= -i(c 2 c 3 - ClC 4 ) 

/oo 
d£ - (pvd^v-x) 

-00 



= -iy / 2\q\Bv(c 2 c 3 - C1C4) 
skB^ 



-i2\q\Bv 1 



(B39) 



and 



/oo 
dy u( a)t (k) 7 °a :E ^ (s) (k) = 2s(A - skP) . (B40) 
-00 

With this we can write down our final expression for the 
transverse pressure for charged particles 



Pi = 



k\B 2 

2tt 2 



/OO 1 
dk z — f + (E k ,T,n) 
-00 ^k 



s=±l n 



|g|z/m(v) 



x/m 2 + 2i/|g|B 
where m(z/) = -\/m 2 + 2z/|g|.B — sk.B. 



• (B41) 
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2. Uncharged Particles 

We now consider the case of uncharged particles. This 
case is different since the transverse momenta of the par- 
ticles are not quantized. Starting from the Dirac-Pauli 
equation (B7) we make an ansatz for the bi-spinor \& of 



the form W(x) — e lk ' x it with u = (ci C2 C3 C4) 1 - . This 
results in the following matrix equation for u 



- kB 





k z 







m + kB 


k+ 


-k z 


k z 


fc_ 


—m + kB 





k + 


-k z 





— m — 



u = Eu , 
(B42) 

where k± = k x ± ik y . Evaluating the determinant of the 
matrix on the left we obtain the energy eigenvalues 



E s = ±^/k 2 + (X-snB) 2 , 



(B43) 



where now we have A = \J m 2 + k\ with k\ = k 2 + k 2 . 
Once again the choice of an overall positive sign corre- 
sponds to particle states and negative sign to anti-particle 
states. We focus on particle states since the result is 
straightforward to extend to anti-particles. 
The resulting positive energy solutions are 




(B44) 



where as before a s = E s — nB + s\ and /3 S = X + sm. The 
overall normalization of the state is fixed in this case by 
requiring that = 2E s S rs . The general quantum 

state for positive energy states can now be constructed 



li Jk 



1 



s=±l ■ 



l2E k 



6 s (k)u( s) (k)e a > a 



(B45) 



where b s (k) is a particle creation operator and J k — 
(2ir)~ 3 f d 3 k. Once again the factor of y/2Ek in the 
denominator above is fixed by the spinor normalization 
used above. 

Following the same general procedures used in the 
charged particle derivation one obtains the following re- 
sult for the energy density 



„ li J k 



s=±l 



The result for the number density is 

n = (N) = J2 f f+(E k ,T,n). 

8=±l Jk 

The result for the parallel pressure is 

P\\=(T ZZ )= E / ff + (E k ,T^). 



(B46) 



(B47) 



(B48) 



And, finally, the result for the transverse pressure P± 

(t xx ) = (ryy) is 



^=E 



s=±l 



1 

k E k 



1 k , m 



2 yjm 2 + k\ 



snBm f + (E k ,T, fi) . 

(B49) 

where m = U m 2 + k 2 ^ — skB. 

In all of the expressions above, we can perform two of 
the three integrations by making the following change of 
variables 



k x = \J \ 2 — m 2 cos 4> , 
k y = \J \ 2 — m 2 sin cj> , 
k z = yjE 2 — (A - skB) 2 

Evaluating the Jacobian one finds 
EX 



(B50) 



d 6 k 



y/E 2 - (A - snBf 



dE dXt 



(B51) 



With this change of variables we obtain the number den- 
sity 



1 r°° 

-3 E / dE E f + (E,T, fi) 

L ~ s=±l Jm-SKB 



/ dX- 



y/E 2 - (A - skB) 2 
1 f°° 

-j E / dE E f + (E,T,[i) 



k + skB arctan 

V \ k 



skB — m\ 7r 
" 2 



(B52) 



where k = \J E 2 — (m — skB) 2 . The energy density is 
given by 



1 f°° 

-j E / dE E 2 f+(E, T, fi) 

L ~ s=±1 Jm-snB 



fE+SKB 

x I dX 



y/E 2 - (A - skB) 2 



? n 2 E / 

J s= 4- \ J rn— ski 



dEE 2 f+(E,T,n) 



k + skB arctan 



/ skB — 



m \ 7r 



V k 



(B53) 
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For the parallel pressure we obtain 



1 f°° 

pE+snB 

x / dXX^E 2 - (A 2 -skB) 2 . 

J m 

1 f°° 

24 ^ E / dEf + (E,T,n) 

^ 47r ~T, Jm-SKB 



-E l 



J tti—skB 

x ^2k(sKB - m)(2m + sk_B) 

/ / skB — m\ 7r 

4k + 6s^i? arctan — 

V V k 



(B54) 



and for the perpendicular pressure we obtain 

Y poo rE+SK,B 

Pi = ^E / dEf+(E,T,n) / 

z7r „ i i Jm—suB Jm 



dX 



s=±l 



X — skB 



y/W - (A 2 - SKB) 2 



(A 2 - to 2 ) - skBA 



-t 



6tt 2 



s=±l 



.—skB 



dE /+ (£, 7» (fc 3 - 3snBmk) . 

(B55) 
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